Homogenous and heterogeneous magnetism in (Zn,Co)0 
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A series of (ZnO) m (CoO) n digital alloys (m = 2, 8, n = 1) and superlattices (m = 80, n = 5, 10) 
grown by atomic layer deposition has been investigated by a range of experimental methods. The 
data provide evidences that the Co interdiffusion in the digital alloy structures is sufficiently effi- 
cient to produce truly random Zni-^Co^O mixed crystals with x up to 40%. Conversely, in the 
superlattice structures the interdiffusion is not strong enough to homogenize the Co content along 
the growth direction results in the formation of (Zn,Co)0 films with spatially modulated Co con- 
centrations. All structures deposited at 160° C show magnetic properties specific to dilute magnetic 
semiconductors with localized spins S = 3/2 coupled by strong but short range antiferromagnetic 
interactions that lead to low temperature spin-glass freezing. 

It is demonstrated that ferromagnetic-like features, visible exclusively in layers grown at 200° C 
and above, are associated with an interfacial mesh of metallic Co granules residing between the 
substrate and the (Zn,Co)0 layer. This explains why the magnitude of ferromagnetic signal is 
virtually independent of the film thickness as well as elucidates the origin of magnetic anisotropy, 
as observed by us and others. Furthermore, in films grown at 200° C and above we observe a 
superparamagnetic contribution. In this case a sizable nanoparticle magnetic moment originates 
from ferromagnetic metallic Co or Co rich nanoparticles dispersed in the bulk of the layer. 

Our conclusions have been derived for layers in which the Co concentration, distribution, and ag- 
gregation have been determined by: secondary-ion mass spectroscopy, electron probe micro-analysis, 
high-resolution transmission electron microscopy with capabilities allowing for chemical analysis; x- 
ray absorption near-edge structure; extended x-ray absorption fine-structure; x-ray photoemission 
spectroscopy, and x-ray circular magnetic dichroism. Macroscopic properties of the layers charac- 
terized by the above techniques have been investigated by superconducting quantum interference 
device magnetometery and microwave dielectric losses allowing to confirm the important role of 
metallic inclusions. 

PACS numbers: 75.50 Pp, 68.55.Nq, 75.75.Cd, 81.15.Gh 



I. INTRODUCTION 



Since the theoretical suggestion by ab initio computa- 
tions that (Zn,Co)0 can be intrinsically ferromagnetic,— 
and the subsequent experimental observation of high- 
temperature ferromagnetism, 2 this compound has 
reached the status of a model system for a broad class 
of dilute magnetic oxides (DMOs) and dilute magnetic 
semiconductors (DMSs), in which a robust ferromag- 
netism is observed despite a minute concentration of 
magnetic impurities. However, despite the apparent 
agreement between experiment and theory for many of 
such materials systems, it was rather soon realized that 
the origin of the abundant high-temperature ferromag- 
netism is far from being understood.^ Indeed, over the 
recent years it has become more and more obvious that 
the understanding of these ferromagnets requires the use 



of advanced nanocharacterization tools in order to asses 
how magnetic impurities are actually incorporated and 
distributed depending on the growth conditions and co- 
dopingj^ On the theoretical front^ it has been argued 
that strong correlation, disorder, and errors in the band 
gap, to mention only few challenges, make questionable 
the direct applicability of standard ab initio methods to 
these compounds. 6 

Our viewi£ (not necessarily shared by all groups) in 
the context of the present work can be formulated as 
follows. 

First, in the case of a random TM distribution and in 
the absence of valence band holes, no ferromagnetism 
is expected above ~10 K^iLIL In fact, in many stud- 
ies, including our own of (Zn,Mn)Q (Refs. [9j-fl3h and 
of (Zn,Co)0 (Refs. M, IH 51 and [HI) only a paramag- 
netic response has been observed, affected actually by 
antiferromagnetic coupling between neighboring spins. 
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This coupling is in fact so strong that the extrapolated 
Neel temperature would be close to 1000 K in wz-CoO. 
Such a robust antiferromagnetism, contradicting some ab 
initio predictions^ has indeed been found in colloidal 
nanocrystals of wz-Zni-^Co^O studied up to x = 1»« 

Second, the abundant observations of high tempera- 
ture ferromagnetism in DMSs and DMOs, if not origi- 
nating from experimental artifacts^ are brought about 
by a highly non-random distribution of transition metal 
(TM) ions, introduced to the sample either purposely 
or via contamination. 4 This non-random distribution is 
driven by a significant contribution of d levels to bonding 
which results in the formation of TM-rich nanocrystals 
characterized by high ordering temperature. 

In the case of (Zn,Co)0, three origins of heteroge- 
neous high temperature ferromagnetism have been con- 
sidered. First, transmission electron microscopy ) 18 > 19 x- 
ray diffraction^ and x-ray magnetic circular dichroism 21 
give an evidence for the presence of metal Co inclusions. 
Also, the presence of superparamagnetic behavior in Co- 
implanted ZnO has also been assigned to Co nanopar- 
ticles, whose size could be changed by annealing^ In 
our recent preliminary work, by employing depth pro- 
filing x-ray photoelectron spectroscopy, we have found 
that metallic Co is located at the interface to the sub- 
strate in films showing ferromagnetic features^ Sec- 
ond, it has been suggested that the actual structure of 
relevant nanoparticles might be more complex, for in- 
stance, they could consist of an intermetallic ferromag- 
netic CoZn compound. 24 Finally, it has been argued 
that uncompensated spins at the surface of antiferromag- 
netic wz-CoO nanocrystals could give rise to spontaneous 
magnetization^ the effect already observed in the case 
of NiO nanoparticlesi2£ 

In this paper, we present results of magnetic stud- 
ies and of microwave dielectric losses measurements- 
sensitive to metallic inclusions — carried out on (Zn,Co)0 
samples characterized by a number of complementary 
element-specific probes. Our films have been grown by 
atomic layer deposition (ALD) according to a procedure 
described in Sec. [TX| The ALD is a self-limiting growth 
process, introduced in the 1970s by Suntola and Antson 
(see Refs. and l28l ) which had initially only limited 
niche applications (such as, e. g., thin film electrolumi- 
nescence devices) but recently has proven itself highly 
successful and of rapidly growing technological impor- 
tance, particularly when the Intel Company started to 
use ALD for deposition of a high-K oxide (Hf02 and re- 
lated oxides) as a gate dielectric in MOSFET transistors. 
This method not only allows superior control of layer 
thickness and perfect coating of surfaces with different 
shapes and morphology, including the 3D ones, but, and 
most importantly for the present studies, allows a sub- 
stantial reduction of the temperature of the deposition 
process. 

When applied to (Zn,Co)0, ALD provides a high qual- 
ity homogenous nanocrystalline structure with grains ex- 
hibiting a various degree of the orientation of the wurtzite 



(wz) c-axis. Furthermore, this technique allows to grow 
superlattice-like structures (ZnO) m (CoO) n with various 
average nominal values of Co concentrations x nom = 
n/(n + m) and periods n. As detailed in Sees. UlllVIl 
the Co concentration, distribution, and aggregation have 
been determined by: secondary-ion mass spectroscopy 
(SIMS); electron probe micro-analysis (EPMA); x-ray 
absorption near-edge structure (XANES); extended x- 
ray absorption fine-structure (EXAFS); x-ray photoelec- 
tron spectroscopy (XPS); high-resolution transmission 
electron microscopy (HR-TEM) with capabilities allow- 
ing for chemical analysis, and x-ray circular magnetic 
dichroism (XMCD). In Sec. IVIII we present results of 
magnetic measurements carried out employing a super- 
conducting quantum interference device (SQUID) mag- 
netometer. Section IVIIII contains results of microwave 
dielectric loss studies. 

This set of experiments allow us to establish the re- 
lation between the growth parameters and the distribu- 
tion of Co. Namely, when the growth process is car- 
ried out at 160°C then independently of other growth 
details, layer thickness, and Co concentration the layers 
exhibit paramagnetic properties, that is in a sense that 
the other, 'legacy', DMSs did: (i) the level structure and 
magnetism of single TM ions can be described (includ- 
ing the c-axis related magnetic anisotropy) by the rele- 
vant group theoretical model; (ii) their paramagnetism 
is weakened by strong antiferromagnetic superexchange 
among the nearest neighbor TM cations with no evi- 
dence for ferromagnetic coupling for any pair distances, 
and (iii) the samples freeze to a spin-glass state on low- 
ering temperature. In particular, in the case of digi- 
tal alloy superlattices (ZnO) m (CoO)„, n = 1, the tem- 
perature and field dependencies of magnetization as a 
function of the Co content x determined from the Curie 
constant show the behavior expected for randomly dis- 
tributed spins coupled by short range antiferromagnetic 
interactions up to x — 40%. This indicates that interdif- 
fusion homogenizes the Co distribution in these digital 
structures. In contrast, long period (ZnO) rrl (CoO) rl su- 
perlattices (n — 5 or 10) show properties which can be 
ascribed as ZnO/(Zn,Co)0 superlattices. We conclude 
that interdifiusion is not strong enough to randomize the 
Co distribution along the growth direction in the long 
period superlattices. 

Furthermore, we have found growth conditions lead- 
ing to a ferromagnetic-like behavior. More specifically, 
our results demonstrate that its presence is associated 
primarily with the growth temperature. When the 
growth temperature is risen to 200° C or above - the 
ferromagnetic-like response appears. It comes about in 
two main guises. 

Firstly, as a robust, nearly temperature independent 
(specified by high Tc) and highly anisotropic response 
which we unambiguously associate with a few nm thin 
metallic Co mesh located at the (Zn,Co)0/substrate in- 
terface. Its existence is widely documented throughout 
the paper and it explains why the magnitude of the 
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ferromagnetic-like signal is virtually independent of the 
film thickness as well as elucidates the origin of magnetic 
anisotropy, as observed by us and others. 29 Furthermore, 
it makes it possible to understand significant deviations 
from the standard superparamagnetic behavior visible in 
our samples as well as in many high-temperature DMSs 
and DMOs^ 

Secondly, anomalous magnetic response is seen in the 
form of a superparamagnetic-like behavior, pointing to 
the presence of nanoparticles magnetized internally up 
to above room temperature. We show by comparing the 
Co concentration from SIMS data with the one resulting 
from magnetic measurements that the relevant nanopar- 
ticles consist rather of ferromagnetic Co (or intermetallic 
compounds) than of uncompensated spins at the surface 
of aniferromagnetic CoO. Furthermore, at least so far, 
our studies do not provide hints for defect- or hydrogen- 
mediated ferromagnetic interactions. Also, our data do 
not confirm the recently suggested^! relation between fer- 
romagnetism and density of grain boundaries in ZnO. 



II. SAMPLES 

A. Growth method 

The (Zn,Co)0 samples have been grown by ALD at 
substrate temperature between 160°C and 300°C em- 
ploying the F-120 Microchemistry reactor and a double 
exchange mechanism of the growth. We use diethylz- 
inc (DEZn) or dimethylzinc (DMZn) as a zinc precur- 
sor, cobalt(II) acetylacetonate (Co(acac)2) as a cobalt 
precursor and deionized water as an oxygen precursor. 
These highly reactive precursors are sequentially intro- 
duced to the growth chamber, so they meet only at the 
surface of the him. This very specific characteristic of the 
ALD technique means that our samples are deposited as 
(ZnO) m (CoO) n periodic structures with various combi- 
nations of m and n values. We grow our layers in either 
a digital alloy fashion (m = 2or8;n = l)or employing 
a superlattice concept (m = 80 and n — 5 or 10). The 
studied films have been deposited on silicon, glass, and 
sapphire substrates, however as we find no qualitative 
differences among them we concentrate on Si-substrate 
based films. 

Apart from above, a wide range of other parameters 
controls the ALD growth process - the length of the ALD 
pulses (the dwell time of the particular precursor in the 
growth chamber), the waiting time before the purging 
process by an inert gas, time of purging (by nitrogen 
here), growth temperature, etc. This enables the prepa- 
ration of films with quite different crystallographic order, 
stoichiometry and Co content. 



B. Investigated samples 

In Table U we display pertinent parameters character- 
izing growth conditions and properties of samples inves- 
tigated in this study. 

We find that it is the growth temperature of 200°C 
and above that plays the decisive role in developing of 
ferromagnetic-like features in our layers. Independently 
of the growth protocol, thickness of the layers, and an 
average Co content ferromagnetic and / or superparamag- 
netic signatures are present there. According to SIMS 
results presented in Fig. Q] for sample F307, a Co accu- 
mulation near the interface to substrate is visible in such 
samples. This is corroborated by TEM, XPS, XMCD 
data discussed in subsequent sections, which point to the 
presence of metal Co mesh at the interface. 

The layers obtained at 160°C do not show even slight 
signs of ferromagnetic or superparamagnetic features, 
they are as good paramagnet down to spin-glass freez- 
ing temperature as the Co-Co antiferromagnetic superex- 
change allows them to be. A quantitative analysis pre- 
sented in Sec. IVIll demonstrates a perfectly random dis- 
tribution of Co cations, with a clear distinction for m — 
80 and n = 5 or 10 supperlattices, where insufficient out- 
diffusion into ZnO slabs results in the ZnO/(Zn,Co)0 su- 
perstructure instead of a uniform mixed alloy. This mod- 
ulated Co distribution is beyond SIMS spacial resolution, 
which shows a uniform Co content along the growth axis 
for sample F72 (Fig. [J). 

Structural characteristics of the (Zn,Co)0 films has 
been assessed by XRD measurements using the X-Pert 
MRD Pro Alphal diffractomctcr (Panalytical) , equipped 
with an incident beam Ge monochromator and a strip 
detector. Diffraction data are collected in a broad 29 
range (20° - 80°), which includes main ZnO diffraction 
peaks corresponding to (10.0), (00.2), (10.1), and (11.0) 
crystallographic orientations. 

All the samples show a homogenous wurzite-type 




Depth ( urn ) 

FIG. 1. (Color online) SIMS Co depth profile for sample F72 
and F307 deposited at 160° C and 200° C, respectively. A Co 
accumulation at the interface to substrate is visible for sample 
F307. 
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TABLE I. List of the samples investigated in this study. We indicate the growth temperature, the ZnO to CoO cycles' ratio, 
thicknesses, wurtzite c-axis arrangement, Co concentrations obtained from SIMS, EPMA, EDX, and SQUID (when applicable), 
and the Hall electron concentration. This is followed by the established layers magnetic character, when PM stand for purely 
paramagnetic layers fSec. lVIIB"]! . SP denotes layers where in addition to PM a sizable superparamagnetic component dominates 
at elevated temperatures (Sec. IVII C[) . and FM indicates layers which additionally show a strong and temperature independent 
ferromagnetic response fSec. IVII Dl) . 



structure with a various degree of directional ordering 
of the wz c-axis. We find no traces of foreign phases, 
cobalt oxides in particular, at least down to the base 
sensitivity of our equipment, estimated to be about 1%. 
Such an exemplary XRD spectrum is shown in Fig. [2] 
Here, a strong (00.2) reflex indicates that the majority 
of the (Zn,Co)0 wz grains are oriented along the growth 
direction, whereas two weak (10.0) and (10.1) reflexes in- 
dicate that a certain, small, fraction of the grains assume 
random orientations. From each spectra like this using 
Schrerrer's formula—^ we calculate an average grain size 
which stays within 5 to 7 nm for most of our layers. 
However, it assumes considerably larger values of nearly 
30 nm in the case of F72 and F73 m — 80 superlattices, 
where thick slabs of ZnO constitute a vast part of the 
layer. This is in line with our previous findings for ALD 
grown (Zn,Co)0 layers^ and indicates a detrimental role 
of Co incorporation onto the crystallographic fidelity of 
(Zn,Co)0. 

On the other hand, the ALD method allows to control 
the preferred orientation of the c-axis with respect to the 
substrate plane i 12 ' 34 As displayed in Table Q] and shown 
in Fig. [3] textured films with dominant orientation of the 
c-axis along the growth direction (e. g., F53 and F176), 
perpendicular to it (e. g., F72), or polycrystalline films 
with randomly oriented c-axis (samples F175 and F179) 
can be prepared. 

There are however many conditions that influence the 
crystallographic orientation assumed by the grains. In 
addition to our previous studies of ZnO, where the purg- 
ing time and growth temperature were identified as the 
decisive factors^ the present studies clearly show that 



an increase of the water pulsing time (the dwell time of 
the oxygen precursor in the ALD growth chamber) ad- 
ditionally favors the perpendicular orientation of grains 
whereas the presence of cobalt counters any grain order 
since most of the highly Co-concentrated samples show 
no preferred orientation. This is clearly a multidimen- 
sional parameter space and more dedicated studies are 
on the way to further refine this issue. However, as the 
c-axis orientation determines the type of the magnetic 
anisotropy exerted by the layers we come back to this 
issue in Sec. IVII B 5l showing. in particular, that it is pos- 
sible to extract a quantitative information on the c-axis 
distribution using low temperature magnetometry. 
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FIG. 2. (Color online) XRD spectrum of layer F176. Indices 
of crystallographic directions corresponding to the diffraction 
maxima are indicated. 
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FIG. 3. (Color online) XRD spectra of four (Zn,Co)0 lay- 
ers showing ordered (F53 and F72) or random (F175 and 
F179) grains' orientation. The indices of crystallographic di- 
rections corresponding to the diffraction maxima as well as 
corresponding orientation of the wurtzite c-axis with respect 
to the sample plane are indicated. 



III. XANES AND EXAFS INVESTIGATIONS 

XANES and EXAFS experiments provide detailed 
and quantitative data on the local structure, there- 
fore from this studies we extracted information about 
atomic configuration around Co ions and about a sub- 
stitutional and/or interstitial position of cobalt in the 
ZnO matrix. We perform comparable studies on a set of 
paramagnetic and ferromagnetic layers, represented be- 
low by F254 (paramagnetic) and F328 (ferromagnetic) 
(Zn,Co)0 films. 



A. Methodology 

XANES and EXAFS measurements at the K edge of 
Co have been performed at DESY-Hasylab (Cemo sta- 
tion) at liquid nitrogen temperature using a 7-element 
silicon fluorescence detector. Hard x-ray photons used 
in the experiment enabled penetration of the whole layer 
thickness. The description of the experimental results is 
performed employing the IFEFFIT data analysis pack- 
age making use of the Athena and Artemis codes^ The 
passive electron reduction factor S02 = 0.8 is estimated 
from fitting the first shell of the data for the sample F254, 
where Co substitutes Zn. 

Experimental results are compared to ab-initio com- 
putations carried out employing the FEFF 8.4 code^ 
in order to determine how Co atoms are located in the 
lattice. A cluster of 10 A radius is first created by 
using the known crystallographic data for possible Co- 
containing structures and then the XANES spectra are 
calculated employing the XANES, SCF (Self-Consistent 
Field) and FMS (Full Multiple Scattering) cards. The 
Hedin-Lundqvist exchange and correlation potential is 
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FIG. 4. (Color online) XANES experimental spectra com- 
pared to results of computations for metallic Co and various 
compounds containing Co. Spectra are shifted vertically for 
clarity. 

adopted. Several structures are considered: cobalt ox- 
ides, ZnCo204, metallic cobalt, as well as (Zn,Co)0, 
where Co substitutes Zn. 



B. Experimental results - XANES and EXAFS 

The experimental XANES spectrum compared to ab 
initio results is presented in Fig. H] It can be noticed that 
for the sample F254 the measured spectrum shows the 
best agreement with the computed results for (Zn,Co)0. 
We conclude, therefore, that in the case of this sample 
Co atoms occupy cation substitutional positions. For 
the F328 sample the spectral features are slightly weaker 
indicating that a part of Co atoms can be in the form of 
other compound. 

In Fig. El the Fourier transforms of the EXAFS spectra 
and their fits are collected for the two studied samples. 
Fits parameters are displayed in Table HU As seen, the 
atom configurations around Co ion significantly differs 
in these two layers. In the case of the sample F254, the 
first two maxima can be well described by coordination 
spheres expected for ZnO, i. e., four oxygen atoms at the 
distance of 1.97(2) A and six zinc atoms at the distances 
of 3.19(1) A and 3.23(1) A, respectively. 

However, in the F328 layer a new peak appears be- 
tween the two maxima corresponding to the ZnO lattice. 
It is a clear indication for the presence of another crys- 
tallographic phase. In this case, in order to describe the 
experimental results, in addition to ZnO, other possible 
Co environments have been considered. It has been found 
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FIG. 5. (Color online) Magnitude of the Fourier transform 
of the EXAFS spectra and fitting results for the samples. 
Spectra are shifted vertically for clarity. 



Parameter 


F254 


F328 


-Rco-o [A] (4 atoms) 


1.97(1) 


1.98(2) 


SSq 


0.005(2) 


0.007(3) 


-Rco-Zni [A] (6 atoms) 


3.19(1) 


3.19(2) 


-Rco-Zn2 [A] (6 atoms) 


3.23(1) 


3.23(2) 


SS Zn 


0.008(1) 


0.009(2) 


y [%] 




27(8) 


-Rco-Co [A] (12 atoms) 




2.49(2) 


SSco 




0.005(3) 


i?-factor 


0.034 


0.021 



TABLE II. The fitting parameters to EXAFS data shown in 
Fig. [5] Here R denotes the distance to the central atom; ss 2 
is EXAFS Debye- Waller factor; y is the relative contribution 
of metallic Co, and ii-factor indicates the quality of the fit. 

that 73% of (Zn,Co)0 and 27% of metallic Co gives the 
best description of the experimental data. 

To conclude, XANES and EXAFS experiments show a 
different Co local structure for paramagnetic and ferro- 
magnetic (Zn,Co)0 samples. In paramagnetic films Co 
ions occupy substitutional Zn positions, whereas in ferro- 
magnetic ones part of Co ions substitute zinc and a part 
of them are present in another crystallographic phase. 
Comparison with theoretical fitting indicates that the 
latter ones are in a metallic form. To establish where 
metallic Co is located in the layer we have employed XPS 
and XPS profiling, discussed in the next section. 



IV. XPS INVESTIGATIONS 

We have performed the XPS studies in order to deter- 
mine the chemical state of cobalt in (Zn,Co)0 films and 
to establish the location of Co within the layer. The XPS 
spectra are recorded on a PHI 5000 VersaProbe™ scan- 
ning ESCA Microprobe using monochromatic A1-K Q ra- 



diation (hv = 1486.6 eV) from an x-ray source operating 
at 100 /im spot size, 25 W and 15 kV. The high-resolution 
XPS spectra are collected with the analyzer pass energy 
of 23.5 eV, the energy step size 0.1 eV and the photo- 
electron take-off angle 45° with respect to the surface 
plane. Shirley background subtraction and peak fitting 
with Gaussian-Lorentzian-shaped profiles are performed 
for the high-resolution XPS spectra analysis. Binding 
energy (BE) scale is referenced to the C Is peak with BE 
= 284.6 eV. 

XPS is a surface sensitive technique in which the prob- 
ing depth depends on the kinetic energy of photoelec- 
trons, but in any case does not exceed few nanometers^ 
In order to get information on the chemical state of cobalt 
within the bulk of investigated films we apply a depth- 
profiling XPS investigation in which XPS spectra have 
been recorded after sequential etching of the films by Ar + 
ions. A few (Zn,Co)0 films, both paramagnetic (PM) 
and ferromagnetic (FM), have been analyzed in this way. 
In Fig. |5] and Fig. [7] we present the representative re- 
sults for two (Zn,Co)0 films: ferromagnetic (FM), (sam- 
ple F328) and paramagnetic (PM), (sample F179). The 
first 15 nm of the etched film was removed using 0.5 kV 
Ar ion etching with the rate of 1.5 nm per minute, and 
then the sputter rate of 15 nm per minute was used (2 
kV Ar ion etching). Both Ar + sputter rates have been 
measured using the Si02/Si reference sample. 

As a result, we obtain elemental depth profiles that are 
presented in Fig.[B][a) and Fig.^b) for FM and PM films, 
respectively. All profiles reveal the elemental distribution 
in bulk of films from 15 nm depth until the (Zn,Co)0/Si 
interface was reached. The relative atomic concentration 
of zinc, oxygen, cobalt, and silicon are evaluated from 
the intensity of XPS peaks associated with the Zn2p, 
Ols, Co2p and Si2s core levels. Cobalt has been de- 
tected across all investigated (Zn,Co)0 films. However, 
its distribution is found to be different for PM and FM 
samples. In the ferromagnetic (Zn,Co)0 films a substan- 
tial enhancement of the Co concentration is found in the 
(Zn,Co)0/Si interface region, as shown in Fig.[B]Ja). The 
cobalt concentration there is estimated to be more than 
three times larger than in the rest of the sample. Such a 
Co-rich interfacial layer is not observed in the paramag- 
netic (Zn,Co)0 films [Fig. HJb)] . 

In order to elucidate the chemical nature of cobalt 
compounds formed in the (Zn,Co)0 films we analyzed 
the high-resolution XPS spectra of the Co2p, which have 
been taken at sequential steps of depth-profiling process- 
ing: at the surface of the (Zn,Co)0 film, after removing 
of the 15 nm of the (Zn,Co)0 layer from the top and 
at the (Zn,Co)0/Si interface. The last two results are 
presented in Fig. [7] 

The original Co2p 3 / 2 spectra are presented together 
with deconvoluted peaks. The Co2p 3 / 2 and C02P2/2 
spin-orbit contributions are separated by 15.0-15.2 eV. 

The main Co2p contributions are accompanied by the 
broad shake-up satellites detected at larger BE. 38 Follow- 
ing deconvolution analysis of the Co2p 3 / 2 XPS spectrum 



7 



three forms of cobalt compounds can be distinguished in 
both samples at a depth of 15 nm (Fig. bottom). 

The strongest component, situated at BE = 780.3 eV, 
is characteristic for cobalt oxide.— Two forms of cobalt 
oxides are stable in air at room temperature: CoO and 
C03O4. They could be identified by the value of the Co2p 
spin-orbit splitting and by intensity of shake-up satellites 
and their positions relative to the main Co2p peaks^r— . 
The BE of the Co2p 3 / 2 peak reported for both cobalt ox- 
ides are very close.— However, the XPS spectra of CoO 
are known to show strong satellite peaks above the main 
2p line and the C03O4 exhibits only weak satellites in this 
region^ The Co2p spectra recorded during Ar + depth- 
profiling of both samples show strong satellites what in- 
dicate CoO to be a main cobalt oxide component. This 
means that the Co atoms that substitute Zn in the ZnO 
matrix, ZnO:Co, contribute to this component. 

Another Co2p 3 / 2 XPS state is located at BE = 
781.8 eV. Its chemical shift is characteristic of cobalt sur- 
rounded by -OH groups. The broad and high-intensive 
satellite peaks (BE at about 786 eV) are also indicative 
for a coexistence of -OH groupsi 38 i 41 ' 42 

The Co2p 3/2 contribution at the BE=777.8 eV can be 
assigned to metallic cobalt^ This observation indicates 
that some Co atoms in the (Zn,Co)0 layers form metallic 
clusters and that these clusters are present both in fer- 
romagnetic and paramagnetic films. Based on magnetic 
data it appears that in samples deposited at low temper- 
atures the concentration and/or the size of Co inclusions 
is too small to be visible in magnetization measurements. 

A significant difference between paramagnetic and fer- 
romagnetic films is observed in the region of (Zn,Co)0/Si 
interface (see Fig. [7J top). In FM film we observe mainly 
metallic cobalt accompanied by very small concentration 
of CoO in the interface region. The absence of shake- 
up satellites, which is a characteristic feature of XPS 
spectra of metallic transition metals and rare earth, indi- 
cate the metallic Co to be a dominant chemical state and 
the cobalt oxide contribution is very small. In the PM 
(Zn,Co)0 film a metallic contribution is accompanied by 
cobalt oxides and Co-OH compounds, but also in this 
case the contribution from metallic cobalt is larger than 
in the volume of the sample. 
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FIG. 6. (Color online). XPS sputter depth-profiles of the a) 
ferromagnetic, b) paramagnetic (Zn,Co)0 films. 



FIG. 7. (Color online) Co2p core level XPS spectra mea- 
sured during Ar + sputter profiling after removing of 15 nm 
of the (Zn,Co)0 film (bottom) and at the (Zn,Co)0/Si inter- 
face region (top). Left: ferromagnetic (sample F328), right: 
paramagnetic (sample F179) (Zn,Co)0 film. Deconvoluted 
spectra indicate different chemical states of Co compounds. 



During the Ar + depth-profiling experiment we have 
also recorded the Ols spectra (not shown here). Three, 
well-separated XPS states can be distinguished. The 
most intensive one at BE ~ 530 eV is attributed to oxy- 
gen bound to zinc. The peak at BE equal to 531.3 eV can 
be related to oxygen bound to Co. The lowest intensity 
peak, located at BE = 532.1 eV, can be assigned to the 
-OH groups. The apparent difference in the shape of the 
Ols peak of FM and PM films is observed in the interface 
region. The Ols XPS spectrum recorded from PM film 
consists two nearly equaled contributions, one related to 
Co-0 bonds, and second to -OH groups. For FM films 
oxygen bound to -OH groups dominates the XPS spec- 
tra. It is an open question what is the role of hydroxyl 
groups in the observed ferromagnetism of (Zn, Co)0. 

In conclusion, the analysis of XPS spectra showed 
that inclusions of metallic cobalt clusters have been de- 
tected in both paramagnetic and ferromagnetic (Zn,Co)0 
samples. However, a different Co-cluster distribution 
has been observed within the (Zn,Co)0/Si interface re- 
gion of both films. The interface in the ferromagnetic 
(Zn,Co)0 film is formed mainly by accumulated metallic 
Co, whereas in the paramagnetic films the metallic Co is 
accompanied by cobalt oxides and hydroxides. 



V. HR-TEM INVESTIGATIONS 

Since both SIMS and XPS studies point to a presence 
of Co enriched regions of our layers we perform detailed 
high resolution transmission electron microscopy (HR- 
TEM) studies in cross section at 300 kV electron beam 
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FIG. 8. High resolution transmission electron microscopy images of a) F179 sample with (Zn,Co)0 wurtzite particle projected 
in the (01.1) direction, and b) F328 sample with Co fee nanoparticle projected in the (Oil) direction. 



energy with the use of a Titan 80-300 Cubed Cs image 
corrected microscope equipped with an energy dispersion 
x-ray (EDX) spectrometer allowing for chemical analysis. 
The cross sectional specimens have been prepared by ion 
milling proceeded with mechanical dimpling. Only digi- 
tal alloy (ZnO) m /CoO layers have been investigated by 
TEM. All the samples reveal a 2 - 3 nanometer thick 
amorphous Si02 layer covering the Si(001) substrate, on 
which the (Zn,Co)0 films are grown. 

The HR-TEM microstructural analysis allows us to 
classify our layers into two main categories. The first one, 
characteristic for layers grown at 160°C, exemplified in 
Fig. [5Ji for layer F179, exhibits a uniform structure with 
oval, 5 to 10 nm in size, monocrystalline grains. These 
values compare very well with the average grain size es- 
tablished from XRD. All these grains show a wurzite 
structure characteristic of ZnO films as can be seen in the 
inset Fig. [8^,. This homogenous polycrystalline structure 
is observed at first 40 - 50 nm from the layer/Si interface. 
In the next part of the film the columnar growth takes 
place and the width of the columns ranges from 25 nm to 
30 nm. No Co-rich volumes are found. According to our 
extended magnetic studies (Sec. IVII[) such layers exhibit 
paramagnetic properties typical for legacy DMS. 

The second category, characteristic for layers grown 
above 160°C, exhibits smaller oval crystallites (diame- 
ter of 3 - 4 nm) at the first 40 - 50 nm of the layer 
and narrower nanocolumns (10 - 22 nm in width) in 
the rest of the layer. However, the basic difference be- 
tween these two types of samples lays in an existence of 
a layer of 3 to 4 nm in diameter Co clusters located at 
the layer/substrate, as exemplified in Fig. [5b for layer 
F328. The close-up on one of such clusters presented in 
the inset to Fig.[S}3 reveals the fee crystallographic struc- 
ture characteristic of metallic cobalt. So, the HR-TEM 
characterization fully confirms conclusion derived from 



XPS studies about the large quantity of metallic cobalt 
present at the interfacial region of such films. It is shown 
later (Sec. IVIID[) that exactly this layer of Co clusters 
that coalesce to a form of metallic mesh spreading all over 
the layer/substrate interface is responsible for a robust, 
nearly temperature independent and highly anisotropic 
ferromagnetic response. 

The point chemical analysis of other layer of this type 
presented in Fig. [9] gives evidence of the existence of Co 
rich volumes also in the bulk of the layer. Samples like 
this one, as presented in Sec. IVII C[ exhibit in an addi- 
tion to the ferromagnetism stemming from the Co-rich 
interface region a clear superparamagnetic-like response 
postulated (see Sec. IVII C|) to originate from about 2 - 
3 nm in diameter hypothetical Co nanocrystallites. We 
immediately note that such objects lay on the very border 
of the detection limit of our crossectional TEM method 
what probably accounts for a lack of their direct visual- 
ization in our HR-TEM images. But, it is equally possi- 
ble that these small Co precipitations assume a coherent 
structure of the wurzite surrounding. 



VI. XAS AND XMCD INVESTIGATIONS 

A. Experimental procedure 

We employ XMCD to verify in the chemical-sensitive 
fashion the origin of the observed ferromagnetic response 
in F328 sample at room temperature. The x-ray ab- 
sorption spectroscopy (XAS) and x-ray magnetic circular 
dichroism (XMCD) measurements have been performed 
at 11011 beamline at MAX-lab synchrotron radiation fa- 
cility in Lund, Sweden42 This beamline is using an ellip- 
tically polarizing undulator source, producing soft x-rays 
of adjustable degree of polarization. It delivers high flux 
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FIG. 9. (Color online) A high angle annular dark field image of the sample F307 with indicated electron beam positions and 
corresponding to them EDX spectra. The cobalt dominates on (Zn,Co)0/Si interface and exhibits high fluctuation inside the 
layer. 



and high brightness circularly polarized x-rays in the en- 
ergy range 0.2 to 1.7 keV, covering the L-edges of the 
late 3d elements. An ultra high vacuum end station has 
been employed in combination with this beamline, allow- 
ing for work in ultra high vacuum (UHV). The samples 
have been mounted in the UHV chamber and character- 
ized at room temperature by means of XAS and XMCD 
after baking the system for about 8 hours at 130°C, to 
reach the base pressure of 5 ■ 10~ 10 mbar. A set of ro- 
tatable coils allows to apply a uniform magnetic field of 
325 Oe, either in the direction perpendicular or almost 
parallel to the film plane, irrespectively of the x-ray angle 
of incidence. The magnetic coils are also employed to de- 
liver magnetic field pulses of 650 Oe, between which the 
XMCD signal can be measured in the state of magnetic 
remanence. The measurements are performed in the total 
electron yield (TEY) mode by measuring the photocur- 
rent of the sample. Using TEY, typically the near surface 
region, from the sample surface down to 6 nm below the 
sample surface is probed with the samples investigated 
here. 



B. XAS and XMCD experimental results 

First, by using the elemental specificity of XAS the 
composition of the (Zn,Co)0 sample is probed for these 
sample positions. Here linearly polarized x-rays are used 
yielding a higher photon flux. The resonant character of 
the excitation makes XAS a more sensitive tool than non- 
resonant XPS or Auger-ES in the laboratory for checking 
the presence of impurities in small concentrations. The 
samples do not contain traces of Ti, V, Mn, Cr, Fe atoms. 
The L edges of these elements fall in the energy range we 
probe. The sample measured "as grown" in the first few 
nanometers contains about 6% of Co ions on average. 
This amount can be calculated by using the tabulated 



atomic cross sections for XAS. The distribution of Co 
is inhomogeneous. The highest amount of Co measured 
on this sample is about 14%, which is close to the value 
obtained from the analysis of the XANES data. Also the 
fine structure of the Co white lines exhibits differences 
versus the "as grown" sample (Fig. ITTJl). 

The Co L-edge for the "as grown" film and around 
the crater shows the typical multiplet shape observed 
for (Zn,Co)0 magnetically diluted samples.— Also the 
O K-edge is similar to the spectra shown in literature for 
(Zn,Co)0. However within the crater both the Co L-edge 
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FIG. 10. (Color online) L-edge x-ray absorption spectra ver- 
sus photon energy for (Zn,Co)0 from the interface region (re- 
gion A, in the crater) and in the "as grown" state (region B 
of the as grown film); the Co L-edge data are taken at normal 
x-ray incidence (a 90° angle between the x-ray propagation 
direction and the film surface plane) in the total electron yield 
mode and normalized to the atomic continuum at high photon 
energies. 
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and O K-edges exhibit differences and indicate both a dif- 
ferent stoichiometry as well as different electronic state 
for the O and Co ion cores probed by XAS. More in de- 
tail, by following the high energy atomic continuum it is 
seen that both the amount of Co and O does decrease 
within the crater. However the inter-peak continuum of 
final states for the Co spectra appears to increase within 
the crater. These final states, which exhibit s symmetry, 
are clearly visible in the spectrum in the inter-peak re- 
gion for metallic Co and can be used as a probe for the 
degree of metal character. The amount of the metallicity 
of the final states appears to increase within the crater, 
indicating that the Co may form there either a contin- 
uous ultrathin film or even small particles. The shape 
of XAS spectra obtained for Co L-edge indicates that 
in the interface region the metallic Co is the dominant 
phase but not the only one. We observe the superposi- 
tion with small amount of Co04£ The difference of the 
electronic state of the Co atoms at the surface versus the 
interface region manifests itself also in the big difference 
in the number of holes of Co atoms in these two regions 
of the sample, illustrated as difference in the area un- 
der L3 line (Fig. ITU)) . The number of holes for Co atoms 
in the surface region is about 7, in the crater is about 
3.5, when for metallic Co this value is 2.8. This result 
confirms the difference in the hybridization of Co with 
neighboring atoms. Inside the crater we not only observe 
a small amount of Co but also a small amount of Zn. 
A small XMCD dichroic signal is seen in the interface 
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FIG. 11. (Color online) XAS (left scale) and XMCD (right 
scale) spectra in the total electron yield mode versus pho- 
ton energy. Measurements are performed at room tempera- 
ture under an applied magnetic field of H — 350 Oe. The 
XAS spectra are obtained with nearly fully circularly polar- 
ized light (a light helicity of 0.85); a 40° x-ray incidence angle 
is employed. The insert shows an enlarged view of the L3 
edge of the XAS and XMCD spectra. 



region (inside the crater) at room temperature, when the 
magnetic field of 350 Oe is applied in the surface plane. 
The dichroic response is much more pronounced at the 
L3 white line indicating that the orbital moment carried 
by the Co atoms at the interface is much stronger than 
for bulk Co. The magnitude of m;/m s obtained by em- 
ploying the magnctooptical "sum rules" for (Zn,Co)0 is 
0.31, whereas for metallic Co typical value is about 0.08. 
This is possible if the Co atoms agglomerate in the form 
of small particles with a metallic Co core . 46 ' 47 Alterna- 
tively, a very rough Co film in the crater with metallic 
patches is also consistent with our measured values. 

Turning now to the overall value of the magnetic mo- 
ment within the crater we note that at 300 K we obtain 
a value of the spin moment of 0.21 /is per Co atom, after 
correcting for the angle of x-ray incidence and the helic- 
ity of the x-rays. This value is of the order of 13% of 
the spin moment for Co. We note in this context that 
a continuous film of only 2-3 atomic layers of Co should 
yield the full magnetic moment at 300 K, under the con- 
ditions of the present experiment. The low value of the 
magnetic moment which we obtain here is a further in- 
dication of stronger finite size effects than in the case of 
a thin film, indicating eventually Co particle formation. 
We conclude that both the high orbital to spin moment 
ratio as well as the low value of the magnetic moment 
after the application of the magneto-optic sum rules in- 
dicate the presence of Co particles or a discontinuous Co 
film within the crater area. 

The value of the magnetic moments determined here 
are distinctively different from those found in ZnO con- 
taining Co particles of about 5 nm in diameter dispersed 
in the bulk of the films. 4 & This difference illustrates the 
fact that the precise size, form and geometric shape of 
inclusions play a crucial role here and determines the 
outcome of an experiment. Conversely, an accumulated 
data base of such studies on well structurally character- 
ized systems may later prove invaluable in the identifica- 
tion of the entities present in the investigated samples. 



VII. SQUID INVESTIGATIONS 

A. Experimental procedure 

Having determined the Co distribution in our 
(Zn,Co)0 layers we turn to investigations of macroscopic 
magnetic properties. In particular, we look for the na- 
ture and magnitude of interactions between dilute Co 
spins as well as for macroscopic manifestations of the 
Co aggregation observed at the nanoscale. According 
to the data presented below, magnetization M(T, H) of 
samples grown at 160°C shows a paramagnetic (PM) be- 
havior, which can be entirely described assuming that Co 
ions occupy only Zn-substitutional positions and that the 
spin-spin interaction is an antiferromagnetic short range 
superexchange, as in canonical II- VI Co-based DMSs. 
On the other hand, in samples grown at higher tempera- 
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tures, in which nanocharacterization reveals metallic Co 
aggregation, two additional contributions to M(H) can 
also be observed: (i) a fast saturating, highly anisotropic 
and temperature independent ferromagnetic (FM) com- 
ponent and (ii) a relatively slowly saturating superpara- 
magnetic (SP) term. The manifestations of PM, FM, and 
SP contributions are discussed in subsequent sections. 

All magnetic results presented here have been obtained 
with SQUID magnetometer (MPMS XL5 of Quantum 
Design) in its base temperature range: from 1.85 K (de- 
noted here for simplicity as 2 K) up to about 390 K fol- 
lowing the procedure described recently^ In order to 
minimize contamination originating from tools and han- 
dling the specimens are rather cleft than cut. Routinely, 
we investigate about 5x5 mm 2 specimens of a roughly 
square shape which are mounted either on a single 20 cm 
long and 1.5 mm wide silicon strip or between two such 
strips for in-plane and perpendicular to plane SQUID 
measurements, respectively. The sample holding strips 
are customary cut from a commercial silicon wafer, pro- 
vide adequate rigidity and stability, and make the mea- 
surements basically artifacts free at the whole accessi- 
ble temperature 2 < T < 400 K and magnetic field H 
(±50 kOe) range. To save on the apparatus time and 
resources (helium), after careful checking for correctness, 
most of the field dependent measurements are performed 
during a one-way field sweeps only, and the missing re- 
turn field sweep and so a full hysteresis loop (if needed 
for presentation) is created by numerical inversion of this 
one-way-sweep data set. When shown, this artificial data 
are marked differently. 

All the films studied here have been deposited on a 
silicon substrate whose magnetic response constitutes 
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FIG. 12. (Color online) Nonlinear part of the magnetic re- 
sponse at room temperature of paramagnetic (Zn,Co)0 films 
on Si substrates. The magnetic field has been swept only once 
for each sample, from about +40 to -30 kOe. The combined 
diamagnetic signal from the substrate and paramagnetic of 
free Co cations has been derived for each sample by the slope 
of its m(H) at large magnetic fields and subtracted from the 
original data. The compensation level for these measurements 
approaches 99.5% already at H = 10 kOe. 



a dominant part of the measured magnetic moment of 
the whole structure. We have found that the substrate 
magnetization consists of the expected diamagnetic com- 
ponent linear in the magnetic field and of an unex- 
pected but rather tiny nonlinear contribution, shown in 
Fig. [TU Whether this sigmoidal dependencies, suggest- 
ing some contamination by ferromagnetic nanoparticles, 
is intrinsic to the Si wafers employed or, less probably, 
results from a sample handling procedure, is unknown 
at present. The shape of this signal is additionally per- 
turbed by a trapped magnetic field in windings of the 
SQUID superconducting magnet M- A typical magnitude 
of this field after returning back to "zero" from 50 kOe, is 
20 Oe. For a usual magnetic response of substrates, this 
field produces a spurious signal of the magnitude compa- 
rable to envisaged values of remanence in DMS films. 

Magnetization values of (Zn,Co)0 films presented in 
the subsequent subsections have been obtained by sub- 
tracting from the total signal the substrate contribution, 
measured as a function of the magnetic field and temper- 
ature, and adjusted for a particular sample by its weight 
and the presence of the trapped field. 



B. Magnetization of Zn-substitutional Co ions 

1. Theoretical modeling 

Following previous extensive studies of various DMSs 
and DMOs ffi 10 ' 49 we exploit magnetic measurements in 
order to evaluate the concentration distribution of Co 
ions as well as the character and magnitude of the ex- 
change interactions between them. In this section we 
summarized theoretical models and quote values of ma- 
terial parameters employed in subsequent sections to 
describe theoretically the field and temperature depen- 
dence of magnetization M(T, H) brought about by Zn- 
substitutional Co ions. 

Properties of a single Zn-substitutional Co 2+ ion (d 7 ) 
in wz-ZnO can be described by the general 5 = 3/2 
Hamiltonian^S: 

H = g\\t-L B H z S z + g^fi B (H x S x + H y S y ) + DS 2 Z1 (1) 

where the z axis of the coordinate system coincides with 
the hexagonal c-axis and H is an external magnetic field. 
The energy level positions are calculated by numerical di- 
agonalization of the 4x4 Hamiltonian matrix (S Z \H\S' Z ) 
with \S Z ) = | - 3/2), | - 1/2), |l/2), |3/2>. Having eigen- 
values, the partition function Z and then magnetization 
M(T, H) the concentration xNq of magnetic ions is ob- 
tained as, 

M(T, H) = k B TxN d In Z/dH, (2) 
This general formula is equivalent to: 

M x<v = fi B g±(S x ,y) 1 M z = Hb9\\(S z ), (3) 
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where 

_ E^i(yjl^l^i) ex p(-^7 fc B r ) (4) 

E- =1 exp(-^/fe B T) 

Here and </?j- are the j-th energy level and the eigen- 
state, respectively. This procedure allows us to calculate 
magnetic response of a Co 2+ ion, which owing to a non- 
zero orbital momentum associated with the spin S — 3/2 
of Co 2+ cations, depends on the angle 8 between the 
magnetic field and the c axis. Electron paramagnetic 
resonance (EPR) ; 50 ' 51 , optical absorption^ 2 - and direct 
magnetometr y 15 ' 53 ' 54 prove the existence of a sizable low 
temperature magnetic anisotropy in Z^^Co^O. 

We describe our data employing the values of pa- 
rameters determined experimentally for Zni^Co^O: the 
cation concentration No = 4.2 x 10 22 cm~ 3 , the Lande 
factors 3|| = 2.238 (H || c), g± = 2.276 (H _L c), and the 
uniaxial single-ion anisotropy energy D = 0.342 meV (~ 
4 K) i 15 ' 53 The positive sign of D results in an 'easy-plane' 
configuration, that is (Zn,Co)0 magnetizes faster with H 
applied perpendicularly to the c-axis of the crystal. 

In order to take into account in the theoretical mod- 
eling of M(H) the presence of short range antiferromag- 
netic spin-spin interactions, we follow the well established 
procedure i 55 ' 56 It consists of replacing x and T in Eq.[5]by 
two temperature dependent fitting parameters, x e R and 
T + TaFj respectively, where Taf > and x c s < x. Typ- 
ically, at helium temperatures, where the nearest neigh- 
bor spins form antiferromagnetically coupled magneti- 
cally inert singlet pairs, x c ff = x(l — x) zi , where z\ = 12 
is the number of the nearest neighbor positions in the 
first coordination sphere in the cation wz-sublatticei££ 

As noted in Sec. Ill B[ our films show a distribution of 
c axis orientations. For polycrystalline samples, with a 
uniform distribution of c axis, we average M(H) over the 
full angle, corresponding to all orientations of the c axis 
with respect to the magnetic field direction. 

In order to describe textured films, in which there is 
a partial ordering of c axis directions, we introduce an 
additional fitting parameter y that describes the fraction 
of the crystal grains having their c-axis perpendicular 
to the sample plane. In terms of this parameter, the 
magnitude of magnetization for the in-plane magnetic 
field assumes the form, 

M(H || plane) = y ■ M(H _L c) + (1 - y) ■ (M(H))g, (5) 

where (M(H))g denotes magnetization average over the 
in-plane angle 9 between the direction of the magnetic 
field and uniformly distributed orientations of the c axis 
in the film plane. 

For the magnetic field perpendicular to the sample 
plane we employ 

M(H ± plane) =y-M(H || c) + (1 -y) -M{H _L c). (6) 

Next, we consider the temperature dependence of mag- 
netic susceptibility. The presence of spin-spin interac- 



tions can be taken into account within the high temper- 
ature expansion^ This procedure leads to the Curie- 
Weiss law which, for a random distribution of localized 
spins, is parameterized by two constants, Co and Go in- 
dependent of the Co concentration x, 

x(T) = xC /(T-xOo), (7) 

where 

©o = -~S(S + l)J2zjJj- (8) 

3 

Here, the summation extends over the subsequent cation 
coordination spheres; Zj is the number of cations in the 
sphere j, and Jj = Jij is the corresponding Co-Co ex- 
change integral in the Hamiltonian fly = —JijSiSj. We 
note that in another convention twice smaller values of 
Jj are considered and then the factor 1 /3 in Eq. [5] is re- 
placed by 2/3. Furthermore, an effective nearest neigh- 
bor exchange energy J c g is sometimes introduced in the 
literature, in terms of which 

9 = ~S{S + 1)Z1 Jeff. (9) 

The values of the Lande factors quoted above point to 
an average value of the Lande factor (g) = 2.25, suitable 
to describe polycrystalline samples. For this magnitude 
of (g) we obtain the Curie constant C — xCo, where 
Co = 0.17 emuK/cm 3 . 

The high temperature expansion is valid as long as 
T 3> |© |, where the Curie- Weiss temperature = xQo. 
In a wider temperature range, x(T) of random antiferro- 
magnets is well described b y 58 ' 59 x{T) = aT~ a , where a 
is a temperature independent constant and a < 1. 

As we show below, this time-honored model of mag- 
netism in II- VI DMSs describes quite satisfactorily dig- 
itally doped samples grown at sufficiently low tempera- 
tures. It also applies to superlattice-like films provided 
that the modulated character of the Co distribution is 
taken into account. 



2. Temperature dependence of magnetic susceptibility 

Figure [T2] presents room temperature dependence of 
(volume) magnetization M(H) on the applied field H for 
selected (Zn,Co)0 layers, where M is obtained from m 
using the whole thickness of the structure to assess the 
volume of the investigated layer. All data show a good 
linearity in H , except the low field region where a resid- 
ual sigmoidal signal, not completely removed by the sub- 
traction procedure described above, mar the otherwise 
paramagnetic M(H). Despite this little technical glitch, 
these data demonstrate that ferromagnetic or superpara- 
magnetic contributions are negligibly small in these films. 
Samples exhibiting such properties are denoted as " PM" 
in Table HI We note that if a sample exhibit such a linear 
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M(H) at room temperature then no traces of ferromag- 
netic component are seen down to the base temperature 
of the magnetometer, T « 2 K. In other words, if a fer- 
romagnetic response is present, it persists to well above 
room temperature independently of a nominal Co con- 
centration. 

We start our quantitative analysis from magnetic sus- 
ceptibility x(T), whose temperature dependence is ob- 
tained from M(H,T) data measured below 10 kOe. The 
values of x(T) determined in this way are plotted in a 
doubly logarithmic scale in Fig. [14] The data indicate 
that xCO T~ a , where a < 1 and its magnitude de- 
creases with shown in the inset. Such a dependence 
is characteristic for random antiferromagnets, that is for 
paramagnetic compounds with a wide spectrum of (anti- 
ferromagnetic) exchange integrals . 58 ' 59 The upper bound 
of this spectrum, i.e., an effective nearest- neighbor ex- 
change integral J e ff can be evaluated by analyzing x" 1 
vs. T in terms of the high temperature expansion, Eq. [7] 

As shown in Fig. [TH the inverse of x(T) points to a 
negative sign of the Curie- Weiss temperature 9, which 
reconfirms the antiferromagnetic character of spin-spin 
interactions. Slopes of l/x(T) vs. T dependencies pro- 
vide the values of Co concentrations x listed in Table Q] 
Furthermore, according to Eqs.[7]and[51 the extrapolated 
values X~ 1 (r = 0) for randomly distributed spins should 
be independent of x and directly proportional to the ex- 
change integrals characterizing spin-spin coupling. And 
the data gathered in the main part of Fig. [T5] instruct 
us that this is the case for the layers grown in the digi- 
tal alloy fashion, indicating that independently of the Co 
content such a growth mode leads to films with randomly 
distributed Co cations over the ZnO host lattice. 
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FIG. 13. (Color online) An example of room temperature 
M{H) for some of the paramagnetic (Zn,Co)0 layers. The 
arrow indicates the direction of the growth of the Co content 
in these layers according to the data accumulated in Tabled] 
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FIG. 14. (Color online) Temperature dependence of magnetic 
susceptibility. The data sets are labeled by the Co concen- 
tration x, established for digital alloy layers according to a 
method described in Sec. IVIIB 31 and by their periods for 
two superlattices. The dashed line indicates the Curie law 
X(T) oc 1/T. Solid and doted lines indicate that \{T) oc T" Q , 
where a < 1, the dependence specific for a random antiferro- 
magnet. The values of a established in this way for the upper 
bunch of curves are plotted vs. x in the inset. The dashed 
black line shows a oc x -1 ^ 3 , and serves as a guide for the eye. 
This trend is not obeyed the samples grown in a superlattice 
fashions, for which the magnitude of \ is low (x ~ 0.7%), 
whereas a (marked as dotted line in the inset) corresponds to 
x w 6%. 



3. Determination of Co concentrations 

Turning now to magnetic assessment of an average Co 
concentration we want to start from a comment concern- 
ing a general reliability of the values of C and 9 de- 
termined from the high-T susceptibility. Our numerous 
experimental and numerical tests showed that whereas 
the slope of % -1 (T) changes substantially in response 
to small variations of the compensating (substrate) data 
that the credibility of 9 value established by this method 
is much greater than that of C. So, in order to estab- 
lish more trustworthy values of x in our layers we re- 
fer back to the room temperature high-field M(H) data 
which are far less sensitive to various experimental arti- 
facts and procedure inaccuracies than the determination 
of the high temperature susceptibility. Therefore we take 
Orj from Fig. [15] and substitute the high-field slope of 
the M(H) in Fig. [13] for \ in Eq. E Obtained this way 
values of x are listed in Table U as our prime magnetic 
estimates of Co concentration in the paramagnetic lay- 
ers. Noteworthy, we find both sets of 'magnetic' x fairly 
equivalent each other (most of the discrepancies are well 
contained within 10% experimental error) and that both 
compare favorably to other methods' estimates (SIMS, 
EDX, and EPMA). The only noticeable discrepancy is 
seen for F254, the most Co-concentrated layer, for which 
9 ~ —200 K and so the high temperature expansion 
(Eq. [7]) ceases to be valid, as the representation of the 
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data via the Curie- Weiss law is valid only when |0| is 
substantially smaller than the temperature of the mea- 
surement. 

We finally conclude this part stating that the reported 
here agreement boosts our confidence in the overall cor- 
rectness of the applied here SQUID-data analysis, in par- 
ticular that other rendered such an experimental proce- 
dure impossible. 54 Furthermore, the data points to the 
random distribution of magnetic ions, which indicates 
that interdiffusion is efficient enough to homogenize the 
Co distribution digitally doped layers. 



4- Long period (ZnO) m (CoO) n superlattices 

It appears to be an altogether different story with the 
long period m = 80, n — 5 and 10 superlattices. As evi- 
denced in Figs . [T4l and [T5l their x(T) stick out completely 
from those of the layers grown in the digital alloy manner. 
In fact, a simple simulation for F72 sample shows that the 
absolute values of x(T) can be well reproduced assum- 
ing 8:1 partition ratio between pure ZnO and Zni-^Co^O 
sublayers with x = 6% in a rectangular ZnO/Zni_ x Coa;0 
superlattice, indicating that interdifussion is not strong 
enough to homogenize the Co content along the growth 
direction in such long period superlattices. 



5. Magnetic anisotropy 

We examine samples for which XRD measurements 
(Fig. [3] in Sec. IIIBI) have revealed the wealth of vari- 
ous c-axis arrangements. These are textured films with 
the dominant orientation of the c-axis either along the 



growth direction (F53) or perpendicular to it (F72) as 
well as polycrystalline samples (with randomized c-axis 
orientations) of a medium (F175) or a high Co content 
(F179), x — 6.4 and 16%, respectively. 

The experimental data shown in Fig. [15] (symbols) 
for parallel and perpendicular directions of the magnetic 
field with respect to the film plane demonstrate that mag- 
netization anisotropy follows the trend expected theoret- 
ically (lines) for the distribution of c-axis revealed by 
the XRD measurements. In particular, in the textured 
films the magnitude of magnetization is larger for the 
magnetic field perpendicular to the prevailing direction 
of the c-axis [Fig. ITBT a.b)]. The description of these data 
has been carried out following the model presented in 
Sec. lVIIB 11 where the parameter y quantifying the frac- 
tion of the crystal grains having their c-axis along the 
growth direction has been introduced. According to the 
fitting procedure, it attains the value of 0.54 in sample 
F53, for which XRD reveals that the c axis is mostly out 
of the film plane, whereas y — 0.10 in the case of sample 
F72, where the c axis is preferably oriented in-plane of the 
film. In contrast, magnetization is isotropic [Fig. ll6f c.d)] 
in the case of samples that are polycrystalline according 
to the XRD results. 

Finally, we want to underline a value of the presented 
here analysis of the low temperature M(H) as an accu- 
rate method for determining the degree of the crystallo- 
graphic order in polycrystals exhibiting a strong magne- 
tocrystalline anisotropy. One has to however keep it in 
mind that experimental handling of very thin layers is 
tricky at any circumstances, and despite a quite strong 
magnetic anisotropy in (Zn,Co)0, a "slight" error in the 
substrate contribution can sizable change the input data 
for fitting and so invalidate the magnetic assessment of 
the c-axis ordering. 
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FIG. 15. (Color online) The inverse of the magnetic sus- 
ceptibility x as a function of temperature for paramagnetic 
(Zn,Co)0 layers. The main panel groups data for lay- 
ers in which CoO layers were introduced in a digital way, 
whereas the inset presents data for the two layers grown in 
a supperlattice-like fashion. Note the substantially different 
values of x _1 (T = 0) for the two types of samples. 



6. Quantifying spin-spin interactions 

The high temperature data on magnetic susceptibility 
x(T) discussed previously (Fig. IT5"j) demonstrates that 
in the films in question Co ions are distributed ran- 
domly. Furthermore, extrapolation of x _1 (T) to zero 
provides magnitudes of Curie- Weiss temperatures for 
particular samples. Now, assuming that the Co-Co cou- 
pling can be characterized by one exchange energy de- 
scribing the interaction with twelve nearest neighbors, 
for the experimental value of x _1 (T = 0) = — Qq/Cq — 
4200 ± 400 cm 3 /emu we obtain 6 = -700 ± 70 K and 
J G ff = 47 ± 5 K. This value compares very favorably 
with 42, 51, and 53 K obtained for single crystalline 
thin films 5 ^ (when we restore the omitted factor "2" in 
their definition of J nn ), single phase bulk crystalls 1 ^ and 
powders^ (after a necessary adjustment of the material 
parameters to the values adopted in this study), respec- 
tively. 

Interestingly, all these values are nearly 40% smaller 
than that established for bulk wz-(Cd,Co)S (Ref. loll ) 
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FIG. 16. (Color online) Magnetization for a series of samples measured at 2, 5, 15 and 300 K (from top to bottom, indicated 
in panel (d)) for two orientations of the magnetic field: in-plane (solid squares) and perpendicular to the sample plane (open 
circles). Lines show modeling of the data by theory outlined in Sec. IVHB II The solid and dashed lines in panels (a,b) are 
calculated for the mixed anisotropy case considering the in-plane and perpendicular magnetic field, respectively. The fraction 
y of grains with the c axis perpendicular the film plane and the effective Co concentration a; ff(T) are the fitting parameters, 
whose values are displayed in the panels (a,b) and in Table Mil respectively. The dotted lines in panels (c,d) represent the 
isotropic (randomized with respect to c-axis) values of magnetization in polycrystalline films. The orange dashed-doted lines 
represent a considerably improved fit with an effective temperature Taf added to the calculations performed at the lowest 
temperature and displayed in Table HTT1 



and twice smaller than those found in bulk zinc-blende 
(Zn,Co)S and (Zn,Co)Se (Ref.H). This implies another 
dependence of | ©o I 011 the bond length d comparing to 
the case of Mn-based II- VI DMS, where tends to de- 
crease monotonically with d i 10 ' 49 i 63 

The quantitative description of M(H, T) provides also 
the values of x e s and Taf , collected in Table IIIII These 
parameters supply also information on spin-spin cou- 
pling. In particular, the magnitudes of x c g determined 
at 2 K are in agreement with the concentrations of Co 
ions having no other Co as the nearest neighbor (see, 
Sec. IVII B 1[) . This means that there is a strong antifer- 
romagnetic coupling between Co ions occupying any of 
the twelve nearest neighbor positions in the cation sub- 
lattice. At the same time, much lower values of Taf (2 K), 
comparing to high temperature Curie- Weiss Q, indicate 
that coupling to next nearest neighbors is relatively weak, 
which points to a rather short range character of the an- 
tiferromagnetic spin-spin interactions. 

Altogether, the findings accumulated so far support 
those ah initio simulations which do not predict the 
presence of ferromagnetic interactions at any distance of 
Co-Co pairs^i^ 



7. Spin-glass freezing 

Since the very early stage of magnetic studies of DMSs, 
a low temperature cusps on x(T) curves were observed 
over a broad range of magnetic cation concentration*^— 
Despite the lack of competing interactions (at least in 
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TABLE III. Temperature dependence of the effective Co con- 
centration x e g, given here as a ratio to its value x at 300 K, 
for the four samples presented in Fig. 1161 This is followed by 
a line giving the statistically expected fraction of single Co 
cations present in the wurtzite lattice at given x, corrected in 
the case of the superlattice sample (F72) by the period ratio 
(Sec. IVII B 4[) . The bottom line lists the values of effective 
temperatures Taf added to the calculations of M(H, T) for 
two samples with the highest x (F175 and F179) at 2 K. 

the compounds where an effective p-type doping is far 
too low to initiate a carrier mediated ferromagnetic ion- 
ion coupling) but on the account of the presence of 
the other two key ingredients, namely positional disor- 
der and spin frustration, these materials exhibit spin- 
glass characteristics driven entirely by antiferromagnetic 
interactions j^Siii And our two paramagnetic layers with 
the highest Co content show a weak low temperature 
feature that closely resembles previous findings in other 
DMS compounds Figured?] collects zero-field cooled 



16 




Temperature ( K ) 

FIG. 17. (Color online) Main figure: a spin-glass-like be- 
havior of two paramagnetic layers with the highest Co con- 
tent. Darker color: F179, I S 16 %; lighter one: F254, 
x 42 %. Zero field cooled (ZFC) and field cooled (FC) 
measurements performed at 1 kOe are indicated by arrows 
and marked by open triangles pointing towards increasing or 
decreasing temperature, respectively. Full circles mark ther- 
moremanent magnetization, TRM (magnified 5 times for bet- 
ter presentation) , which enables us to establish the spin-glass 
freezing temperature, Tf, most reliably. Inset shows Tj de- 
pendence on the transition metal content: im„ in (Zn,Mn)0 
(open diamonds from Ref. [72j and open square from Ref. [73f ) 
and x in (Zn,Co)S (full hexagons from Ref. [fU and red bullets 
in (Zn,Co)0 (this study). 

(ZFC), field-cooled (FC) and thermoremanence (TRM) 
magnetization for these layers. A typical for a transition 
to a glassy state bifurcation on ZFC and FC curves is 
seen. However a strong paramagnetic background masks 
the effect of freezing (characteristic cusp on the ZFC) 
considerably, therefore the TRM measurement serves to 
accurately establish the freezing temperature Tf for these 
layers. 

An important question arises now whether the pre- 
sented here history dependent effects and a very weak 
magnetic hysteresis developing at T < Tf (not shown) 
are indeed related to the spin-glass-like freezing, or they 
are a manifestation of the dynamical blocking of some 
small ferromagnetic clusters (Co-metal droplets or some 
Co-rich spinel precipitates) that together give rise to 
supcrparamagnetic-like behavior. These doubts are fu- 
eled by our earlier observation that a contribution from 
metallic Co is present in the XPS spectra both of the 
PM layer F179, one of these large x layers which show 
the freezing. In order to resolve this issue one should 
resort to studies of the long time scale dynamics of the 
non-reversible part of the signal— However, the presence 
of a large PM background renders such an experiment 
problematic. 

In an attempt, we have found a lack of any noticeable 
time dependence of ZFC magnetization up to 2-10 4 s after 
the x = 16% sample was zero-field-cooled to 5 K, that is 
below its Tf (not shown) , what is in strike difference with 



fast relaxations observed in ZFC state for superparamag- 
netic (SP) layer F309, presented in Sec. I VII CI however 
on the account of a very small accuracy, we regard this 
finding as rather weak evidence. 

Therefore, we arc going to argue for the spin-glass 
freezing on the account of the minute magnitude of M, 
corresponding solely to the PM response of these two 
samples. We note, that if sizable and numerous Co- 
rich metallic (ferromagnetic) inclusions were present they 
would dominate the magnetic response at weak magnetic 
fields (at any temperature) considerably increasing the 
magnitude of M, see e.g., Ref. p75h . Indeed, as it is exem- 
plified in Sec. I VII C[ the FC magnetization of this afore- 
mentioned SP layer F309 is found be nearly 500 times 
larger at low temperatures than in the presented here 
layers, despite ten times weaker magnetic field used dur- 
ing the measurement of the SP sample. Therefore, pre- 
dominantly on the account on the size of the irreversible 
effects, we rule out their superparamagnetic origin and 
adopt the spin-glass explanation. 

Finally, we find a good correspondence between the 
observed here freezing temperatures with those reported 
for similar systems. As indicated in the inset to Fig. [17] 
found here Tf follow the already established exponential 
dependency T f ~ x a with a = 2.1±0.4 for bulk (II, Co) VI 
compounds^ 4 - and a ~ 2 for the other legacy DMSs . 69 i 76 
Interestingly, Tf data from Ref. [72| do not joint none of 
these trends indicating that an unidentified factor must 
have contributed to the observed there freezing-like be- 
havior. 



C. Bulk superparamagnetic contribution 

We find that the total amplitude of M(T,H) can, in 
general, be decomposed into three main components: (i) 
the already discussed paramagnetic (PM) term brought 
about by randomly distributed Zn-substitutional Co 
cations; (ii) a relatively slowly saturating in the mag- 
netic field superparamagnetic (SP) component, and (hi) 
a fast saturating, highly anisotropic and temperature in- 
dependent up to 300 K ferromagnetic (FM) contribution. 

The SP component is present in all our samples grown 
at 200°C and above. Additionally, in such samples a FM 
contribution may show up. We postulate now, and elab- 
orate it later, that it is few nm thick mesh of metallic 
cobalt present on the layer/substrate interface which is 
responsible for the ferromagnetic component. The exis- 
tence of such a film has already been documented in the 
previous sections devoted to general characterization of 
our layers (SIMS, TEM, XPS, and XMCD). We also re- 
gard this interface FM to be completely temperature in- 
dependent, as finite-size effects are not supposed to show 
up for Co film thickness above ~ 1.5 nm. 77 

Conversely, the PM and SP components are assumed 
to stem from the bulk of the layer. So, we allow all 
the Co atoms which build up the volume of the layer to 
either take random cation sites, giving rise to the PM 
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response, or to aggregate into nano-crystallites of a large 
magnetic 'supermoment' to, responsible for the SP part 
of M(T,H). 

Dotted lines in Fig. [TBI indicate a result of the modeling 
of the experimental data by a straight sum of these three 
magnetic components. We calculate the PM response fol- 
lowing the method described in the previous section. In 
particular, we use isotropic paramagnetic M(T, H), allow 
for the temperature dependence of x e g, and take Tap = 
for simplicity. The magnitude of this signal is parame- 
terized by the concentration of randomly distributed Co 
cations, x G g(300 K). 

As described in the next section, the FM component 
has been derived from the room temperature M(H) of 
the layer F328 and scaled to this sample according to 
their areas, so there is no free parameters here. The fit- 
ting itself is performed only for H > 5 kOe, as we do not 
have enough knowledge on the processes which determine 
coercivity and remanence and hence the magnitude of M 
at weak magnetic fields. 

A priori, as mentioned in the Introduction, spins con- 
tributing to the SP signal might originate from two dis- 
tinct kinds of systems: (i) nanoparticles of metallic Co, 
CoZn, or related intermetallic ferromagnetic compounds; 
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FIG. 19. (Color online) Blocking phenomenon in a super- 
paramagnetic sample (F309). Main part: temperature de- 
pendence of the zero field cooled (ZFC), field cooled (FC) 
and thermoremanent (TRM) magnetization. Arrows indicate 
positions of the blocking temperature (7b) and the maximum 
blocking temperature (T^ ax ). Insert: time evolution of the 
magnetization in these three magnetic states measured at 5 K 
(below T B ). 
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FIG. 18. (Color online) A representative (volume) magneti- 
zation for a sample grown at 200° C (F307). Solid squares 
are obtained with the magnetic field applied in the sample 
plain. For clarity the data obtained in perpendicular config- 
uration are shown only for 300 K (open circles). Dotted lines 
indicate results of the modeling which assumes a presence of 
three independent magnetic contributions: a paramagnetic, a 
superparamagnetic, and a ferromagnetic one, as detailed in 
the text. The dashed line represents the magnitude and room 
temperature curvature of the superparamagnetic component, 
whereas the two solid lines represent a derived from sample 
F328 the anisotropic ferromagnetic component used in the 
modeling (see Sec. IVII Dfl . 



(ii) uncompensated spins residing at the surface of anti- 
ferromagnetic wz-CoO nanocrystals, as observed in the 
case of nanoparticles of NiO^ In the latter case, involv- 
ing magnetically inert Co atoms in wz-CoO nanocrystals, 
the total Co density should be larger than the one seen 
magnetically. 

The SP contribution to M(T, H) is approximated by 
the Langevin function mNspL(x), where x = mH/kp,T 
and to = n[iQ . Guided by XPS and TEM data, as a 
first and successful guess we assume that nanoparticles 
constitute of Co metal. In such a case, /ico = 1-7 [ip, 
and the number n of Co atoms giving rise to the mag- 
netic moment to and the density of nanoparticles TVsp are 
bound together by the requirement that the total density 
of the Co atoms, A^a;eff(300 K) + nNgp is equal to the Co 
concentration found by SIMS in the volume of this layer, 
x SIMS (F307) = 8%. Therefore, out of three quantities, n, 
Nsp, and a; o ff(300 K), we are left with two temperature 
and field independent adjustable parameters. 

The fitting results presented in Fig. [T5] demonstrates 
that 65% of Co present in the bulk of the layer assumes 
ZnO lattice cation sites indicating that this is still the 
preferred site for Co atoms. The remaining 35% of Co 
precipitates into nanoparticles containing (on average) 
n ~ 750 Co atoms (approx. 2.5 nm across) if they in- 
deed assume the form of bulk Co. The magnitude of 
this superparamagnetic component at 300 K is presented 
in Fig. [TH] as a dashed line. As seen, this is the super- 
paramagnetic contribution that is responsible for the pro- 
nounced high temperature curvature and slope, whereas 
the paramagnetic component governs the low temper- 
ature magnitude of M at high magnetic fields. The 
interface-related ferromagnetic component accounts for 
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the quickly magnetizing and anisotropic component. Fi- 
nally, we note that obtained this way such a good repro- 
duction of the experimental data strongly indicates that 
all these three components are to the first order indepen- 
dent one from another. 

For further studies of the SP contribution, particu- 
larly to examine blocking properties, we select a sam- 
ple with no FM component. The relevant findings are 
presented in Fig. [191 where a clear and quite sharp max- 
imum of ZFC magnetization and a bifurcation between 
ZFC and FC magnetizations are seen at about Tq = 7 
and = 22 K, respectively. These blocking tempera- 

tures allow to assess an average and maximum size of sta- 
tistically relevant nanocrystals. If surface anisotropyZS is 
unimportant, we can take magnetocrystalline anisotropy 
density K of bulk cobalt, K — 5 x 10 6 erg/cm 3 . Employ- 
ing the widely accepted formulae for DC type SQUID 
magnetometry, Tb — KV/25kB, where V is the volume 
of the nanoparticle, we obtain the average and the statis- 
tically relevant maximum size of Co nanocrystals as 2.4 
and 3.3 nm, respectively. Both values correspond pre- 
cisely to our previous estimate of 2.5 nm, obtained from 
Co counting in the different layer. 

Finally, we confirm the SP nature of the magnetic sig- 
nal in question by examining relaxation of ZFC, FC, and 
TR magnetizations below Tb, as shown in the inset to 
Fig. [TpJ The lack of noticeable relaxation of the FC signal 
and relatively fast relaxations of both ZFC and TR mag- 
netizations, upwards towards the FC value and down- 
wards towards zero, respectively, agree with the general 
expectations for the time evolution of these three type of 
magnetization in a blocked state of superparamagnets. 



D. Interfacial ferromagnetism 



On the account of the already presented XPS, HR- 
TEM and XMCD studies we unambiguously assign the 
swiftly magnetizing and anisotropic signal to the evi- 
denced there few nm thin layer of metallic Co mesh and 
we will show that indeed it shows many characteristics 
corresponding to bulk cobalt. 

In Fig. [2UJ we plot a 'sheet' magnetization Ma, the 
magnetic moment per unit area, obtained at 300 K for 
layer F328, the thinnest layer of all studied layers and 
the main subject of the aforementioned characterization 
efforts. The dashed line marks in the figure the expected 
saturation level, and these are actually the points taken 
from this line plus the original data taken from — 10 < 
H < 10 kOe field range (for both orientations of H) 
which are used in the modeling of the M(T, H) presented 
as the solid lines in Fig. [T5] Interestingly, similarly to 
layer F307, we can completely remove this ferromagnetic 
component from experimental data in some other layers 
(in both configurations of H) by a simple subtraction 
of these established here prototype M(H), providing we 
scale the data according to the layers' areas first. 



However, the TEM cross-section images of the F328 in- 
terface area (see Fig. [5b and Fig. ^ indicate some kind of 
a granular form of this interfacial Co film. Moreover, our 
XMCD studies (see Se dVI|) of the interfacial part of the 
layer also point to a granual form of the Co layer at the in- 
terface. So, the important question arises whether these 
Co granules are (magnetically) independent or whether 
they coalesce towards a uniform Co layer. In order to dis- 
tinguish between these two possibilities we perform zero 
field cooled (ZFC) and field cooled (FC) measurements 
at weak magnetic field. A presence of a well developed 
maximum on the ZFC curve associated with a lack of 
a corresponding one on the FC magnetization would in- 
dicate the temperature dependent (superparamagnetic) 
blocking phenomenon (as already evidenced for a super- 
paramagnetic layer F309 in Fig. [TP)) , and so it would 
indicate a loose distribution of cobalt granules on the in- 
terface. But, as indicated in the inset to Fig. [501 we nn d 
no real maximum what rules out the individual (super- 
foments scenario in favor for the coalescing picture, and 
so pointing strongly to a percolating mesh formed from 
closely packed and interconnected Co granules. 

Another support for the Co mesh comes from the no- 
tion that when we recalculate Ma onto M (volume), we 
obtain M ~ 500 emu/cm 3 for a reasonable value of the 
interface film thickness, dc a = 4 nm. This is only a third 
of a bulk cobalt saturation magnetization what indicates 
that indeed Co granules occupy only a (small) fraction 
of the interfacial region. We also note that the estab- 
lished granular form of the Co interfacial film does not 
contradict our initial assumption about the lack of finite- 
size effects. In this point we share the view of Ney and 
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FIG. 20. (Color online) Magnetization of layer F328. Main 
part: magnetic moment per unit area, Ma, at 300 K. Solid 
squares are obtained with the magnetic field applied in the 
sample plain, open circles mark the data obtained in perpen- 
dicular configuration. The dashed horizontal line indicates 
the estimated saturation level of the ferromagnetic part of 
the data which we use to construct a prototype data sets for 
this magnetic contribution. Inset presents low field in-plane 
zero field cooled (ZFC) and field cooled (FC) magnetization. 
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coworkers 1 ^ that Co granules which are 3-5 nm big should 
not show any pronounced reduction of the magnetization 
between helium and room temperature. 

Finally, we note that it is the flatness of this Co film 
(mesh) that is responsible for the universality of the easy- 
plane magnetic anisotropy. However, the experimentally 
observed anisotropy field, H^ p ~ 5 kOe is about 3 to 
4 times weaker than the expected for a perfectly flat 
and uniform Co film, H^ 1 ~ 18.5 kOe. We take this 
anisotropy weakening as yet another manifestation of the 
patchy character of the Co film, and treat the correspon- 
dence of the magnitude of H^ p to the magnitude of Co 
magnetocrystalline anisotropy field as purely coinciden- 
tal, predominately due to the lack of crystallographic or- 
dering of the Co granules seen in TEM images. 

VIII. MICROWAVE CONDUCTIVITY 
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FIG. 21. (Color online) Spheroidal model of a two-layer sam- 
ple. The electric field component of the microwave electro- 
magnetic radiation is parallel to the longest axis of the sample. 
The ei is the film dielectric permittivity, £2 is the substrate 
dielectric permittivity. 



The presence of metal Co inclusions visualized by XPS, 
as described above, while invisible in standard conduc- 
tivity measurements, could lead to detectable dielectric 
losses. In particular, the microwave cavity perturbation 
technique, 79 ' 80 that is the observation of modifications of 
the cavity resonance by the presence of a sample, may 
serve to evaluate the magnitude of microwave conductiv- 
ity ctac °f our films and to find out if it correlates with 
their magnetism . 



A. Theory of microwave cavity perturbation 
technique 

We use microwave cavity perturbation technique, first 
developed to measure the complex dielectric permittiv- 
ity of small spherical samples placed in a microwave 
cavity 79 ' 80 This method was then adopted to the case 
of ellipsoidal samples^ and subsequently generalized to 
cover a wide range of conductivity magnitudes, from the 
quasi-static case up to the skin-depth regime . 82 ' 83 In most 
cases the cavity perturbation method has been applied 
to uniform samples only. One of a few exceptions is the 
consideration of a layered structure consisting of a thin 
conducting film deposited on a dielectric substrate^ 

In order to apply this method to thin films of 
(Zn,Co)0, we have developed a more rigorous and, simul- 
taneously, more general theoretical approach. We con- 
sider spheroidal samples consisting of two layers (the film 
and the substrate), as shown schematically in Fig. I2T1 
We adopt the quasi-static approximation 84 ", i. e. we as- 
sume the effective wavelength of the microwave radiation 
within the layer to be much greater than its transversal 
dimension. Moreover, we assume that a relative devi- 
ation from the cavity resonance frequency ujq is small, 
however the perturbed microwave electric fields E\ and 
Ei within particular layers may differ considerably from 
the value Eq for the empty cavity. The fundamental per- 
turbation formula, describing the modification of cavity 



resonance by the presence of the sample, can be then 
concisely written as, 

*-iA/2 = ai[^)-l] + °2[6 2 H-l] (1Q) 
' 1 + ni[ei(w) - 1] + na[e 2 (w) - 1] ' V ' 

where 5 = (ujq — uj)/lu and A = 1/Q — 1/Qo denote 
a relative shift from the values of ujq and the quality 
factor Q brought about by the sample; the indices 1 
and 2 refer to the film and the substrate, respectively; 
on = Vi\E 0s \ 2 /(2 J dv\E \ 2 ), where V$ are the sample- 
related volumes while the integration extends over the 
whole cavity; Eq s is the value of Eq at the sample posi- 
tion; ei(u) = e'i{oj) — ie"(w) are complex dielectric per- 
mittivities; and rii are depolarization factors^ It can be 
easily seen that Eq. (flO|) is a direct generalization of a 
similar formula derived earlier for a single layer M- 

Assuming that the substrate parameters are known, 
one can determine the film permittivity ei(w) from ex- 
perimental values of S and A. For a strongly conduct- 
ing film the expression describing e'i(u>) is numerically 
unstable, i. e., we cannot determine e'i(ui) reliably if 
e"(uj) 3> e'i(w). However, Eq. (flU)) can be effectively 
solved for e"(oj), leading to the expression which is anal- 
ogous to that derived 8 ^ for a single layer sample and con- 
nects directly the cavity bandwidth with the layer con- 
ductivity. For flat samples and highly resistive substrates 
the influence of substrate properties is small and can be 
neglected to a good approximation. 

The above exact analytical solution for layered me- 
dia has been derived for spheroidal samples. In prac- 
tice, however, one deals often with samples in the form 
of discs or rectangular prisms and the depolarization fac- 
tors for sample shapes other than ellipsoidal are known 
only approximately. Furthermore, the theory is devel- 
oped for symmetric structures fFig. !21[) in which the film 
is placed in the symmetry plane of the spheroid, while 
real films are typically deposited close to the substrate 
surface. Therefore, replacing a real sample by an equiva- 
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lent model ellipsoid in data analysis could be a source of 
a systematic error. Fortunately, both numerical analysis, 
and test measurements showed^ that for the film pa- 
rameters distinctly different from those of the substrate 
that error is small and measurement results only weakly 
depend on the film position. 



Results of microwave studies 
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Our measurements are carried out at room temper- 
ature employing a cylindrical cavity of the diameter 
d = 49 mm and height h — 35.5 mm, which operates in 
the TE\yi mode of resonant frequency of ~ 9.2 GHz. The 
sample is placed in the maximum of the electric field in 
the plane perpendicular to the cavity axis at the distance 
I = h/A from the cavity bottom. The resonance curve 
(microwave power reflected from the cavity measured at 
frequencies near and at the resonance) is sampled for 
400 frequency values for the cavity without and with the 
studied sample. The fitting to a standard Lorentzian 
function allows to determine both the frequency shift and 
the bandwidth change. The filling factors ol\ and «2 are 
calculated as the ratio of the film and substrate volume 
to an effective volume of the cavity, dependent on its di- 
mensions and the distribution of the electromagnetic field 
for a given resonant mode. It is assumed that for the fiat 
samples under consideration depolarization factors rij are 
small but nonzero, making then possible to measure also 
highly conducting samples. 

Measurements of microwave AC conductivity have 
been performed for a series of (Zn,Co)0 films. Because 
this method can sensitively detect even small highly con- 
ducting objects we search for correlation between ctac 
and layer's magnetic characteristics, since the latter de- 
pend strongly on the presence and type of Co-rich inclu- 
sions, as it has been elaborated in the previous section. 
Importantly, this is a contactless method, what allows to 
investigate also films which are highly resistive according 
to DC measurements. 

We start the description of the results by noting that 
no evident correlation between the AC electric conduc- 
tivity and the DC one is observed for our (Zn,Co)0 films. 
In particular ctac does not correlate with the free electron 
concentration deduced from the Hall effect, as it is indi- 
cated in Fig. Neither it correlates with Co content 
represented in the same figure by its EDX values, as this 
method provides us with an averaged Co content in the 
whole layer, including the Co-enriched interfacial film. 
More importantly, and to some extent surprisingly, we 
see no correlation with the magnetic properties of the lay- 
ers, however ctac appears to assume the largest values in 
samples exhibiting the interfacial ferromagnetism, or to 
be precise here, in layers having the film of interconnected 
metallic cobalt granules on the layer/substrate interface 
(Sec. IVIID[) . and so providing an effective conducting 
medium within the layer. Three such layers (F268, F328 
and F338) are exemplified in Fig. [521 when we plot sheet 
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FIG. 22. (Color online) Sheet AC conductivity (<tac times 
layer thickness) versus EDX Co content. Labels indicate 
layer's approximate electron concentration (in cm -3 , see Ta- 
ble U for more details). 

AC conductivity (ctac times layer thickness) versus the 
EDX Co content. We note that expressed this way ctac 
changes for these samples only by a factor of 2.5 indi- 
cating an important, if not dominant, role of metallic 
inclusions at the interface in both the determination of 
the magnitude of ctac an d in the observed discrepancy 
between DC and AC measurements. Finally, we note 
that the observed such a large disagreement between DC 
and AC conductivity of these layers confirms our previ- 
ous assessments about a discontinuous character of the 
interfacial Co layer. 

On the other hand layer F309 exhibits a similar mag- 
nitude of the sheet AC conductivity despite the fact it 
shows only the superparamagnetic response fSec. lVH~Cj) . 
In this case the sizeable conductivity stems solely from 
the large volume of the film (about 20 times larger than 
of F328) indicating that indeed the method is capable of 
addressing the tiny and separated metallic (conductive) 
objects present in the investigated specimen. However, 
in this case also, the correct assignment of the experi- 
mental finding would not be possible without the joint 
characterization effort of all the methods employed in 
this study. 

IX. CONCLUSIONS 

A series of (ZnO) m (CoO)„ digital alloys (m = 2, 8, n = 
1) and superlattices (m = 80, n — 5, 10) grown by atomic 
layer deposition have been investigated by a range of ex- 
perimental methods. The data provide evidences that 
the Co interdiffusion in the digital alloy structures is suf- 
ficiently efficient to produce truly random Zni_ K Co^O 
mixed crystals with x up to 40%. Conversely, in the 
superlattice structures the interdiffusion is not strong 
enough to homogenize the Co content along the growth 
direction results in the formation of (Zn,Co)0 films with 
spatially modulated Co concentrations. High resolution 
SQUID measurements have demonstrated that all struc- 
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tures deposited at 160°C exhibit magnetic properties spe- 
cific to dilute magnetic semiconductors with localized 
spins 5 = 3/2 randomly distributed over cation sites, and 
mutually coupled by strong, but short range, antiferro- 
magnetic interactions. The determined magnitude of the 
exchange energy describing coupling between the nearest 
neighbor Co pairs corroborates quantitatively its previ- 
ous determination 5 ^ and qualitatively those first princi- 
ples studies, in which the spin-spin coupling was found 
to be antiferromagnetic— ~— The presence of positional 
disorder within the spin system coupled by antiferro- 
magnetic interactions lead to low temperature spin-glass 
freezing, not observed earlier in (Zn,Co)0. According to 
our findings, characteristics of these freezing are similar 
to those reported for other II- VI DMSs. 

Apart from a possible influence of some other growth 
parameters, it is the growth temperature in excess of 
160°C that is the decisive factor leading to the occurrence 
of ferromagnetic-like features in our layers. Whereas our 
previous investigations 23 indicated a correlation between 
accumulation of Co metal at the interface and the ob- 
served room temperature ferromagnetic response, in the 
present work we have investigated the character of this 
accumulation and verified that indeed the FM response 
comes from the interface region where Co is accumulated. 
The SIMS, TEM, and XMCD investigations have indi- 
cated that the inclusions are in the form of somewhat 
oblate (few monolayer thick) Co particles or a flat mesh 
connected magnetically. The latter is in line with the re- 
sults of magnetic measurements showing clearly that the 
ferromagnetic easy axis is in-plane. 

The Co mesh at the interface may suggests that the in- 
terface acts as a Co sink which stops to operate once the 
mesh is built. This would explain why the magnitude of 
the ferromagnetic signal is virtually independent of the 
thickness of (Zn,Co)0 film. Finally, the assignment of 
the ferromagnetic features to the interfacial layer makes 
it possible to understand significant deviations from the 
classical superparamagnetic behavior. Such properties 
are visible in many high-temperature magnetically doped 
semiconductors and oxides^ However, the Co layer is 
not continuous according to transport investigations, as 
we observed the noticeable difference between the DC 
and AC (microwave) conductivity that indicates the pres- 
ence of highly conductive areas within our (Zn,Co)0 films 
but no DC metallic-like conductivity. 



In addition to the ferromagnetic signal, particularly 
well visible in thin samples, in all layers grown at high 
temperature a superparamagnetic response has been ob- 
served. Guided by XPS data, this magnetization compo- 
nent has been assigned to metallic Co-rich nanoparticlcs 
disperse over the film volume. This view has been further 
supported by comparing the magnitude of the total mag- 
netic signal to the Co concentration provided by SIMS. 
The important question is whether the presence of Co in- 
clusions results automatically in superparamagnetic-like 
features. According to the XPS results this is not nec- 
essarily the case, as some concentration of metal Co in- 
clusions have been detected in all samples, also in those 
grown at low temperature and showing only a paramag- 
netic response. In contrast, in samples grown at higher 
temperature, presumably because of more efficient aggre- 
gation of sizable Co-rich nanoparticles, the superparam- 
agnetic contribution is significant. 

In general terms, our results presented here as well as 
in parallel work on TM-doped nitrides j&Si reemphasize 
the necessity of using a range of nanocharacterization 
tools in order to assess the lateral and vertical distribu- 
tions of TM ions in the studied films. A non-random dis- 
tribution of TM atoms and a non-random distribution of 
TM-rich aggregates underline pertinent features of these 
systems, particularly their magnetic properties. These 
distributions depend sensitively on growth parameters, 
which make it possible to obtain a variety of different 
systems, ranging from DMSs with a random distribution 
of magnetic cations to novel nanocomposites. 
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